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Abstract. In this article we consider the strong approximation of sto- 
chastic difFcrcntial equations driven by Levy processes or general semi- 
martingales. The main ingredients of our method is the perturbation 
of the SDE and the Taylor expansion of the resulting parameterized 
curve. We apply this method to develop strong approximation schemes 
for LIBOR market models. In particular, we derive fast and precise algo- 
rithms for the valuation of derivatives in LIBOR models which are more 
tractable than the simulation of the full SDE. A numerical example for 
the Levy LIBOR model illustrates our method. 

1. Introduction 

The main aim of this paper is to develop a general method for the strong 
approximation of stochastic differential equations (SDEs) and to apply it 
to the valuation of options in LIBOR models. The method is based on the 
perturbation of the initial SDE by a real parameter e and then on the Taylor 
expansion of the resulting parameterized curve around e = 0. Thus, we 
follow the line of thought of Siopacha and Teichmann (2007) and extend 
their results from continuous to general semimartingales. The motivation 
for this work comes from LIBOR market models; in particular, we consider 
the Levy LIBOR model as the basic paradigm for the development of this 
method. 

The LIBOR market model has become a standard model for the pric- 
ing of interest rate derivatives in recent years. The main advantage of the 
LIBOR model in comparison to other approaches is that the evolution of 
discretely compounded, market-observable forward rates is modeled directly 
and not deduced from the evolution of unobservable factors. Moreover, the 
log-normal LIBOR model is consistent with the market practice of pricing 
caps according to Black's formula (cf. Black 1976). However, despite its ap- 
parent popularity, the LIBOR market model has certain well-known pitfalls. 

On the one hand, the log-normal LIBOR model is driven by a Brownian 
motion, hence it cannot be calibrated adequately to the observed market 
data. An interest rate model is typically calibrated to the implied volatility 
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surface from the cap market and the correlation structure of at-the-money 
swaptions. Several extensions of the LIBOR model have been proposed in the 
literature using jump-diffusions, Levy processes or general semimartingales 
as the driving motion (cf. e.g. Glasserman and Kou 2003, Eberlein and Ozkan 
2005, Jamshidian 1999), or incorporating stochastic volatility effects (cf. e.g. 
Andersen and Brothcrton-Ratcliffe 2005). 

On the other hand, the dynamics of LIBOR rates arc not tractable un- 
der every forward measure due to the random terms that enter the dy- 
namics of LIBOR rates during the construction of the model. In particular, 
when the driving process has continuous paths the dynamics of LIBOR rates 
are tractable under their corresponding forward measure, but they are not 
tractable under any other forward measure. When the driving process is 
a general semimartingale, then the dynamics of LIBOR rates are not even 
tractable under their very own forward measure. Consequently: 

(1) if the driving process is a continuous semimartingale cap lets can 
be priced in closed form, but not swaptions or other multi-LIBOR 
derivatives; 

(2) if is a general semimartingale, then even caplets cannot be priced 
in closed form. 

The standard remedy to this problem is the so-called "frozen drift" approx- 
imation, where one replaces the random terms in the dynamics of LIBOR 
rates by their deterministic initial values; it was first proposed by Brace 
et al. (1997) for the pricing of swaptions and has been used by several au- 
thors ever since. Brace et al. (2001), Dun et al. (2001) and Schlogl (2002) 
argue that freezing the drift is justified, since the deviation from the original 
equation is small in several measures. 

The frozen drift approximation is by far the simplest approximation pos- 
sible, however it is well-known that it does not yield acceptable results, es- 
pecially for exotic derivatives and longer horizons. Therefore, several other 
approximations have been developed in the literature; in one line of research 
Daniluk and G^tarek (2005) and Kurbanmuradov et al. (2002) are looking 
for the best lognormal approximation of the forward LIBOR dynamics; cf. 
also Schoenmakers (2005). Other authors have been using linear interpola- 
tions and predictor-corrector Monte Carlo methods to get a more accurate 
approximation of the drift term (cf. e.g. Pelsser et al. 2005, Hunter et al. 
2001 and Glasserman and Zhao 2000). We refer the reader to Joshi and 
Stacey (2008) for a detailed overview of that literature, and for some new 
approximation schemes and numerical experiments. 

Although most of this literature focuses on the lognormal LIBOR market 
model, Glasserman and Merener (2003b, 2003a) have developed approxima- 
tion schemes for the pricing of caps and swaptions in jump-diffusion LIBOR 
market models. 

In this article we develop a general method for the approximation of the 
random terms that enter into the drift of LIBOR models. In particular, by 
perturbing the SDE for the LIBOR rates and applying Taylor's theorem we 
develop a generic approximation scheme; we concentrate here on the first 
order Taylor expansion, although higher order expansions can be considered 
in the same framework. At the same time, the frozen drift approximation 
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can be embedded in this method as the zero-order Taylor expansion, thus 
offering a theoretical justification for this approximation. The method we 
develop yields more accurate results than the frozen drift approximation, 
while being computationally simpler than the simulation of the full SDE for 
the LIBOR rates. 

The method we develop is universal and can be applied to any LIBOR 
model driven by a general scmimartingale. However, we focus on the Levy 
LIBOR model as a characteristic example of a LIBOR model driven by a 
general semimartingale. 

The article is structured as follows: in section 2 we review time-inhomoge- 
neous Levy process, and in section 3 we revisit the Levy LIBOR model. 
In section 4 we derive the dynamics of LIBOR rates under the terminal 
martingale measure and express them as a Levy-driven SDE. In section 5 
we develop the strong Taylor approximation method and apply it to the 
Levy LIBOR model. Finally, section 6 contains a numerical illustration. 

2. TiME-INHOMOGENEOUS LEVY PROCESSES 

Let {Cl,J^,F,'P) be a complete stochastic basis, where = Tt^ and the 

filtration F = (.^t)t6[o,T«] satisfies the usual conditions; we assume that T* G 
M^o is a finite time horizon. The driving process H. = {Ht)Q<t<T» is a time- 
inhomogeneous Levy process, or a process with independent increments and 
absolutely continuous characteristics. That is, H is an adapted, cadlag, real- 
valued stochastic process with independent increments, starting from zero, 
where the law of Ht, t G [0,T*], is described by the characteristic function 

E [e'"^*] = exp I y [ibsu - + J (e^'''^ - 1 - i?xx)F,(dx)] ds j ; (2.1) 

here 6t G M, ct G M^o and Ft is a Levy measure, i.e. satisfies Ft{{0}) = 
and /jjj(l A \x\'^)Ft{dx) < oo, for all t G [0,T*]. In addition, the process H 
satisfies Assumptions (AC) and (EM) given below. 

Assumption (AC). The triplets {bt,ct,Ft) satisfy 

J (^\bt\ + ct + J{1A \x\^)Ftidx)^dt < oo. (2.2) 

R 

Assumption (EM). There exist constants M,e > such that for every 
u G [-(1 + e)M, (1 + e)M] =: M 

J J e"^Ft(dx)di < OO. (2.3) 

{|a;|>l} 

Moreover, without loss of generality, we assume that J^^^^^ij e^^Ft{dx) < oo 
for all t G [0, T*] and u G M. 

These assumptions render the process H = (i?t)o<t<T« a special semi- 
martingale, therefore it has the canonical decomposition (cf. Jacod and 
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Shiryaev 2003, II.2.38, and Eberlein et al. 2005) 

H = J bsds + J ^sdWs + J J xin^ -iy){ds,dx), (2.4) 

OR 

where fi^ is the random measure of jumps of the process H, v is the P- 
compensator of , and W = (Wt)o<t<T« is a P-standard Brownian motion. 
The triplet of predictable characteristics of H with respect to the measure 
P, T(ii'|P) = {B,C, v), is 

B = J bsds, C = J Csds, i^{[0,-]xA) = J J Fs{dx)ds, (2.5) 

OA 

where A G the triplet {b, c, F) represents the local characteristics of H. 

In addition, the triplet of predictable characteristics {B, C, v) determines the 
distribution of H, as the Levy-Khint chine formula (2.1) obviously dictates. 

We denote by Kg the cumulant generating function associated to the in- 
finitely divisible distribution with Levy triplet {bs,Cs,Fs), i.e. for z G M, 

se[0,T] 

Ks{z) := bsz + + / (e"" - 1 - zx)Fs{dx). (2.6) 

R 

Using Assumption (EM) we can extend Kg to the complex domain C, for 
z G C with "Rz G M, and the characteristic function of Ht can be written as 

t 

E [e^^-f^*] =^^vi^j Ks{iu)d^ . (2.7) 



If is a Levy process, i.e. time-homogeneous, then (5^, c^, F^) and thus also 
Kg do not depend on s. In that equals the cumulant (log-moment) 

generating function of Hi. 



3. The Levy LIBOR model 

The Levy LIBOR model has been developed by Eberlein and Ozkan 
(2005), following the seminal articles of Sandmann et al. (1995) and Brace 
et al. (1997) on LIBOR market models driven by Brownian motion; see 
also Glasserman and Kou (2003) and Jamshidian (1999) for LIBOR mod- 
els driven by jump processes. The Levy LIBOR model is a m,arket m,odel 
where the forward LIBOR rate is modeled directly, and is driven by a time- 
inhomogeneous Levy process. 

Let = To < Ti < • • • < T/v < T]v+i = T^^ denote a discrete tenor 
structure where 5i = Tj+i — Ti, i G {0. 1. ... , A^}. Consider a complete sto- 
chastic basis (r2, JT, F, Pt^^) and a timc-inhomogeneous Levy process H = 
{Ht)o<t<Tt satisfying Assumptions (AC) and (EM). The process H has pre- 
dictable characteristics (0, C, u'^*) or local characteristics (0, c, F'^*), and its 
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canonical decomposition is 

H = J ^sdWj*+ J J x{n" -i^'^*){ds,dx), (3.1) 

OR 

where W'^* is a P^^^ -standard Brownian motion, /i^ is the random measure 
associated with the jumps of H and v^* is the Pt, -compensator of jj,^ . We 
further assume that the following conditions are in force. 

(LRl): For any maturity Tj there exists a bounded, continuous, deter- 
ministic function A(-, Tj) : [0, Tj] — >■ M, which represents the volatility 
of the forward LIBOR rate process L{-,Ti). Moreover, 

N 

^\X{s,Ti)\<M, 

i=l 

for all s G [0,r*], where M is the constant from Assumption (EM), 
and A(s, Tj) = for all s > Ti. 
(LR2): The initial term structure B{0,Ti), 1 < 2 < + 1, is strictly 
positive and strictly decreasing. Consequently, the initial term struc- 
ture of forward LIBOR rates is given, for 1 < i < A^, by 

The construction starts by postulating that the dynamics of the for- 
ward LIBOR rate with the longest maturity L(-,Tjv) is driven by the time- 
inhomogeneous Levy process H and evolve as a martingale under the ter- 
minal forward measure Pt, • Then, the dynamics of the LIBOR rates for 

the preceding maturities arc constructed by backward induction; they are 
driven by the same process H and evolve as martingales under their associ- 
ated forward measures. 

Let us denote by Pt^+i the forward measure associated to the settlement 
date Ti+i, i e {0,. . ., N}. The dynamics of the forward LIBOR rate L{-,Ti), 
for an arbitrary Tj, is given by 

L{t,Ti) = L{0,Ti)exp b^{s,Ti)ds + J X{s,Ti)dHj'+'^ , (3.3) 

where H^^+^ is a special semimartingale with canonical decomposition 

t t 
Ht'+' = J ^sdwl'+' + J J x{n" -u'^'+^){ds,dx). (3.4) 



Here W'^^+'^ is a PTj+i -standard Brownian motion and zz-^'+i is the PTj+i" 
compensator of . The dynamics of an arbitrary LIBOR rate again evolves 
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as a martingale under its corresponding forward measure; therefore, we spec- 
ify the drift term of the forward LIBOR process L(-, Tj) as 

b'^{s,Ti) = -^X\s,Ti)cs 

-J {e^^''^'^'' - 1 - X{s,Ti)x)Fj'+'{dx). (3.5) 

The forward measure Pt^^^, which is defined on {'^,J',{J't)Q<t<Ti+i), is 
related to the terminal forward measure Py^ via 

(3.6) 

The Pj^^^-Brownian motion W'^*+'^ is related to the Py^-Brownian motion 



via 



t 

Wf'+^ = Wl''+^ - I a(s, Ti.^)^/^As = ... 



't - '"t -J a(s,rj+i)vL.s 



- / f E (3.7) 

n \l=i+l / 



where 



The Pt^+i -compensator of /U^, is related to the Pt, -compensator of 

fi^ via 

i/^'+i(ds,da;) = p{s,x,Ti+i)u'^'+^{ds,dx) = ... 

= (YI P{s,x,Ti)\ u'^*{ds,dx), (3.9) 
\/=i+l / 

where 

Remark 3.1. Notice that the process iJ^'+i, driving the forward LIBOR 
rate L{-,Ti), and H = H^* have the same martingale parts and differ only 
in the finite variation part (drift). An application of Girsanov's theorem for 
semimartingales yields that the PTi+i -finite variation part of H is 

JcsYl (^{s,Ti)ds + J Jx(y[ P{s,x,Ti) - l\ u^'{ds,dx). 

'='+1 M \l=i+l / 

Remark 3.2. The process H = H^* driving the most distant LIBOR rate 
L{-,Tn) is - by assumption - a time-inhomogeneous Levy process. However, 
this is not the case for any of the processes iJ-^'+i driving the remaining 
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LIBOR rates, because the random terms i'^SiL(t- Ti) c'l^ter into the compen- 
sators z/^'+i during the construction; see equations (3.9) and (3.10) 



4. Dynamics under the Terminal Measure 

In this section we derive the stochastic differential equation (SDE) that 
the dynamics of an arbitrary LIBOR rate satisfy under the terminal measure 
Ft, ; this will be the starting point for the perturbation of the dynamics and 
the Taylor expansion in the next section. Of course, we could consider the 
SDE as the defining point for the model, as is often the case in stochastic 
volatility LIBOR models, cf. e.g. Andersen and Brotherton-RatclifFe (2005). 

Starting with the dynamics of the LIBOR rate L(-, Tj) under the forward 
martingale measure P^^-^^, and using the connection between the forward 
and terminal martingale measures (cf. eqs. (3.7)-(3.10) and Remark 3.1), 
we have that the dynamics of the LIBOR rate L{-,Ti) under the terminal 
measure are given by 



L{t, Ti) = L(0, TO exp I J b{s, Ti)ds + J A(s, T^)dHs j , (4.1) 



where H = {Ht)Q<t<T, is the PT^^-time-inhomogeneous Levy process driving 
the LIBOR rates, cf. (3.1); the drift term b{-,Ti) has the form 

N 



b{s,Ti) = -\x\s,T,)cs-CsX{s,T,) -i^^f_^A(.,TO 

- / /'(e^(«'^^)--l) H p{s,x,Ti)-X{s,Ti)x]Fj'*{dx), (4.2) 

H V l=i+l ) 



where (3{s,x,Ti) is given by (3.10). Note that the drift term of (4.1) is 
random, therefore we are dealing with a general semimartingale, and not 
with a Levy process. Of course, L{-,Ti) is not a Pr_^ -martingale, unless 
j = 1 (where we use the conventions that Yl^=i ~ ^-^d n?=i = !)• 

Let us denote by Z the semimartingale that drives the LIBOR rate L(-, Tj) 
under the terminal measure, i.e. 



Z:= J b{s,Ti)ds + J X{s,T,)^sdW, 





+ / f X{s,Ti)x{ix" -u^*){ds,dx). (4.3) 
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The triplet of predictable characteristics of Z = {Zt)o<t<Tt. under Pr, , 
T(Z|FtJ = (B,C,u), is described by 



B 



I b{s,Ti)di 
Jo 



C= / X'{s,T,)csds (4.4) 

^0 

lAix) *v = \a (A(s, Ti)x) *v^% ^(K\{0}) . 

The assertion follows from the canonical decomposition of a semimartingale 
and the triplet of characteristics of the stochastic integral process; see, for 
example, Proposition 1.3 in Papapantoleon (2007). 

Now, using Lemma 2.6 in Kallsen and Shiryaev (2002), we know that we 
can express the natural exponential of the semimartingale Z as the stochastic 
exponential of another suitable semimartingale i.e. e'^ = ^"(2'), where Z 
has dynamics 

Z = Z+^(Z") + (e^-l-x)*/i^ 



J j \{s,Ti)x{ii^ -v^'-){ds,dx) 

R 

+ / /(e^(^'^^)^-l-A(s,ri)a;)M''(ds,dx) 





+ 





= j b{s,Ti)ds+^J X^{s,Ti)Csds + J X{s,Ti)dHs 



+ J J{e^^''^i> -1- X{s,Ti)x)n^{ds,dx). (4.5) 

R 

Therefore, under the terminal measure Pt, the dynamics of an arbitrary 
LIBOR rate L(-,Tj) are described by the SDE 

dL{t,Ti) = L{t-,Ti)dZt, (4.6) 

where Z is a semimartingale with decomposition 

Z = J b{s,Ti)ds + J X{s,Ti)dH, 



+ J J (e^(^'^^)- - 1 - X{s, Ti)x)ti''{ds, dx), (4.7) 
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and the (random) drift term, comparing (4.2) and (4.5), is 

" S,L{s-,T,) 



M.,T.) = -c,A(.,T.)E_j^^;^A(..T,) 

TV \ 

(^gA(.,T,). _ 1^ JJ f3{s,x,Ti)-X{s,Ti)x\Fj*{dx). (4.8) 



=i+i 

The triplet of predictable characteristics of the semimartingale Z under 
Ft, is T(Z|PtJ = {B^,C^,i^^), where 

B^= f {bis,Ti) + K!:iXis,Ti)))ds 
Jo 

= f X\s,Ti)csds (4.9) 
Jo 

l^{x) * z/^ = U(e^(^'^*)^ - 1) * Ae B{R\{0}). 

Here /t'" denotes the cumulant associated to the jumps of the process H, i.e. 

f>(n) = / (c"^ - 1 - ux)F^*{dx). (4.10) 
Jr 

Note that the martingale part - or equivalently the diffusion and jump char- 
acteristics - of the semimartingale Z corresponds to a Levy process. 

5. Strong Taylor approximation and applications 

The aim of this section is to strongly approximate the stochastic differen- 
tial equations for the dynamics of LIBOR rates under the terminal measure. 
This pathwisc approximation is based on the strong Taylor approximation 
of the random terms in the drift 6(-, Tj) of the semimartingale Z driving the 
LIBOR rates L{-,Ti), cf. equations (4.6)-(4.8). The essence of the strong 
Taylor approximation is the perturbation of the initial SDE and a classical 
Taylor expansion with usual conditions for convergence (cf. Definition 5.1). 

We introduce a real parameter e and will approximate the terms 

L{t-,Ti) (5.1) 

which cause the drift term to be random, by their first-order strong Taylor 
approximation, cf. Lemma 5.4. Note that the map x i— > appearing in 

the drift, is globally Lipschitz with Lipschitz constant 5* = max; Si. 

The following definition of the strong Taylor approximation is taken by 
Siopacha (2006); see also Siopacha and Teichmann (2007). Consider a smooth 
curve e ^-^ We, where e G M and G L'^{n;'R). 

Definition 5.1. A strong Taylor approximation of order n > is a (trun- 
cated) power series 

i=0 

such that 

E[|W,-T"(iy,)|] = o(e"), (5.3) 

holds true as e ^ 0. 
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Then, for Lipschitz functions / : M ^ M^o with Lipschitz constant k we 
get the following error estimate: 

nifiWe) - f{T^{W,m < kn\We - T-(M/,)|] = koie^). (5.4) 

Remark 5.2. It is important to point out that motivated by the idea of the 
Taylor series we perform an expansion around e = and the estimate (5.4) 
is valid. However, for the pathwise approximation of LIBOR rates we are 
interested in the region e 1, and hope that the expansion yields adequate 
results; for e = we would simply recover the "frozen drift" approximation. 
Numerical experiments show that this approach indeed yields better results 
than the "frozen drift" approximation; cf. section 6. 

We assume that the following condition is in force, which guarantees that 
L{t,Ti) e L2(0), for alH G {0, . . . ,iV} and all t G [0,T^]. 

Assumption (L). The volatility functions satisfy |A(s,rj)| < ^ for all 
s G [0,Ti] and ah i£{0,...,N}. 

Let us introduce the process X'^(Ti) = (X^(Tj))o<t<T, with initial val- 
ues X^{Ti) = L{0,Ti) for ah i G {0, . . . ,iV} and an~e~G M. The dynamics 
of X^{Ti) is described by perturbing the SDE of the LIBOR rates by the 
perturbation parameter e: 

dX^{Ti) = eXUTi) {b'{t, T,)dt + X{t, Ti)dHt 

+ J (e^(*'^^)^ - 1 - A(t, r,)x)/x^(dt, dx)) , (5.5) 



where. 



^^^l + 6iXUTi) 

_ f ^e^(t,T,)x _ 1^ Yl P^{t,x,Ti)-X{t,Ti)x\Fl'{dx), (5.6) 
R V l=i+i ) 

and 

Remark 5.3. In the sequel we will use the notation T as shorthand for T^. 

Lemma 5.4. The first-order strong Taylor approximation oj the random, 
variable X^{Ti) is given by: 

T{X^{Ti)) = L(0,r,) + e^J^^^X^iTi), (5.8) 
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where the first variation process of X^(Ti), denoted byY(Ti) := ^\f:=QX'^ (Ti) , 
is a time-inhomogeneous Levy process with local characteristics 



= L(0, T,) Ti) + 4"(A(s, Ti))^ 

cf = (^L{0,Ti)X{s,Ti)ycs (5.9) 

I lA{x)FY'{dx) = I u(L(0,T,)(e^(^'^^)- - l))F^*(dx), A e B{R). 

Proof. By definition, the first-order strong Taylor approximation is given by 
the truncated power series 

T{Xim))=X?iT,) + e^J^^^Xtm). (5.10) 

Since the curves e i-^ X^{Ti) are smooth, and X^(Ti) £ L^(r2) by Assumption 
(L) , we get that strong Taylor approximations of arbitrary order can always 
be obtained, cf. Kriegl and Michor (1997, Chapter 1). In particular, for the 
first-order expansion we have that 

JE[\X^iTi)-T{Xtm))\]=oie). (5.11) 

The zero-order term of the Taylor expansion trivially is 

X^{Ti) = X^{Ti) for all t, since dX^{Ti) = 0. 

Of course, the initial values of the perturbed SDE coincide with the initial 
values of the un-perturbed SDE, hence X^{Ti) = L(0,Tj). 

The first variation process of X^{Ti) with respect to e is derived by dif- 
ferentiating (5.5); hence, the dynamics is 

d(^—\^^^X^{Ti)') = XUTi)\e=o{b'{t,Ti)\,=odt + X{t,Ti)dHt 

+ J (e^(*.^^)^-l-A(t,ri)x)/i^(di,da;)) 

= X^{Ti) (6°(i, Ti)dt + X{t, Ti)dHt 

+ J [e^it,T,)x _ ^ _ ^(^^ ri)x)/x^(dt, dx)) . (5.12) 

We can immediately notice that in the drift term b^{t, Ti) of the first vari- 
ation process, the random terms Xf{Ti) are replaced by their deterministic 
initial values X^{Ti) = L{0,Ti). 

Let us denote by Y{Ti) the first variation process of X^(Tj), i.e. Y{Ti) = 
^\^=QX'^{Ti). The solution of the linear SDE (5.12) describing the dynamics 
of the first variation process yields 

t t 

Yt{Ti) = L{0,T,)(^J b\s,Ti)ds + J X{s,Ti)dHs 



t 

+ J J (e^(^'^^)^-l-A(s,ri)x)/x^(ds,dx)^. (5.13) 
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Since the drift term is deterministic, if is a time-inhomogeneous Levy pro- 
cess and is the random measure of jumps of a time-inhomogeneous Levy 
process, we can conclude that Y is itself a time-inhomogeneous Levy process. 
The local characteristics of Y are described by (5.9). □ 

Concludingly, we have the following strong Taylor approximation scheme 
for the random terms (5.1) entering the drift: 

Xt{Ti) = L{0,Ti) + eYt{Ti) + o(e). (5.14) 

Remark 5.5. A consequence of this approximation scheme is that we can 
interpret the "frozen drift" approximation - first proposed by Brace ct al. 
(1997) and then used by several other authors ~ as the zero- order Taylor 
approximation, i.e. Xf{Ti) ^ L(0,Tj). 

Remark 5.6. Let us elaborate on the random terms of the first variation 
process Y{Ti). Define the process 

U = [ X{s,Ti)dH„ 
Jo 

and consider the exponential transform of U, i.e. the process U such that 
U := £og(e^). Applying Lemma 2.6 in Kallsen and Shiryaev (2002) we know 
that U and U are related via 

[/ = C7+^(C7^C7^) + (e^-l-a:)*/x^. (5.15) 

We immediately observe that up to the multiplicative constant L(0,Tj), 
the martingale parts of U and of the first variation process Y(Ti) (5.13) 
are indistinguishable. Now, by the definition of the stochastic logarithm we 
know that U is the solution of the following stochastic differential equation: 

dU = ^d(e^) ^U= [ ^d(e^). (5.16) 

e^- J e^- 

The latter observation is useful for the discretization and simulation of the 
first variation process. 

Next, we parameterize the Levy LIBOR model in terms of the perturba- 
tion parameter e; that is, the stochastic differential equation describing the 
evolution of the LIBOR rates takes the form 

dL'{t, Ti) = L'{t-, Ti) (b'it, Ti}dt + Xit, Ti)dHt 

+ j (e^(*.^*)^-l-A(i,rOx)M''(di,dx)), (5.17) 

R 

where h^{t, Ti) and /3^(t, x, Ti) are given by (5.6) and (5.7) respectively, hence 
they involve the process X'^(Tj). Moreover, let L'^(0,rj) = L(0,rj) for all 
values of e G R and all LIBOR rates, i.e. for alH G {0, . . . , N}. Notice that 
for e = 1 the processes X'^iTi), L^{-,Ti) and L{-,Ti) are indistinguishable. 

Now, we provide a way to pathwise approximate L^{-,Ti) by means of the 
strong Taylor approximation. 
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Proposition 5.7. Assume that LIBOR rates are modeled according to the 
Levy LIBOR model and let e G M. Then, the stochastic differential equation 
for L^{-,Ti), given by (5.17), (5.6) and (5.7), can be strongly approximated 
as € [0 by 



dL'{t, Ti) = L'{t-,Ti) (b'{t, Ti)dt + X{t, Ti)dHt 

+ j (e^(*'^^)^-l-A(t,ri)x)/i^(dt,dx)), (5.18) 



where 



N 



b'{t,Ti) = -ct\{t,Ti) J2 J 



-Kt,Ti) 



N \ 

A(t,r,)x_i^ ]-J p^^t,x,Ti)-\{t,Ti)x\Fl^{dx), (5.19) 

l=i+l ) 



and 



d^t X T,) = '^''^(^*-(^^))+ (^xm). _ 1 w 1 



(5.20) 



Proof. The first step is to interchange X^(Tj) with X^(Tj)-|-, which yields no 
difference to the dynamics of L''^-, Ti) in (5.17), (5.6) and (5.7), since Xf(Ti) 
is a positive random variable. Then, using Taylor's expansion, we know that 
P{t,Ti) L^{t,Ti) P-a.s as e i 0. 

The error in the approximation stems from the drift term, since the mar- 
tingale parts of L'^{-,Ti) and L'^{t,Ti) are identical. We have the following 
estimate for the error term: 



iogL^(t,r,)-iogL^(i,r,)| 
t 

< J\b'is,Ti)-b'is,Ti)\di 



< j \cs\{s,Ti) [oc'{s,Ti)-a'{s,Ti)) 

\ r ( ^ ^ \ 

+ y ly (e^(^'^^)- - 1) n /^'(^'^^^o- n Hs.^.T{)\Ff{dx) 

Q ]g \/=i+l J=i+1 / 



ds. 
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Now, for the first term we can estimate using the Lipschitz property of 



X 



l+Six 



and Assumption (L): 



N / 



5iXl_{Ti) 5{Y{XUTi)) 



+ 



+ 5iXl_{Ti) 1 + 5{T{XI_{T{))_ 



ds 



< j c,|A(.,r,)|<5* \xUm-T^{XUTi))^\\Ks,Ti)\A. 



/=i+l 



Similarly for the second term, using again the Lipschitz property, Assump- 
tions (EM) and (L), and in addition that 11^=1 (^A; + 1) — 1 + X^fc=i o^k for 
small ajfc, we get 

* „ / N N \ 

j I j (e^(^''^^)- - l) ns,x,Ti)- n ^\s,x,Ti)\fJ*{Ax)\As 

Q ]g \i=i+l l=i+l / 

* N 



t 



l=i+l 

N 



< I k{M){N - i - 1)5* ~T{Xl_{Ti))_ 



l=i+l 

Therefore we can conclude that 

t 



ds. 



N 

I log L^(i,T,)- log L^(t,T,) I < / Q,5* ^ \xl_(Ti) -T{Xl_(Ti)) ^\ds, 
where Qg '■= ^si^, M) is a constant depending on s, N and M. 



□ 



Naturally, we can solve the stochastic differential equation (5.18) describ- 
ing the dynamics of the approximate LIBOR rate L'^{-,Ti). Let Z'^ denote 
the semimartingale 



Z' = J b'{t,Ti)dt + J \{t,Ti)dHt 





+ / / (eMW)--l-A(t,ri)x)/i^(dt,dx); 
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T 


0.5Y 


lY 


1.5Y 


2Y 


2.5Y 


B{0,T) 


0.9833630 


0.9647388 


0.9435826 


0.9228903 


0.9006922 


T 


3Y 


3.5Y 


4Y 


4.5 Y 


5Y 


5(0, T) 


0.8790279 


0.8568412 


0.8352144 


0.8133497 


0.7920573 



Table 6.1. Euro zero coupon bond prices on February 19, 2002. 



then we get that 

dP{t,Ti) = L'{t-,Ti)dZt 



L\t,Ti) = L\d,Ti)£{Z% 

= L^(0, Ti) exp (Zl - + (log(l + x) - x) * iif) 



t t 



L^(0, Ti) exp J b\s, Ti)ds + J \{s, Ti)dHs^ , (5.21) 



where 







b\s, Ti) = b'{s, Ti) - ^A2(s, ri)c„ (5.22) 



and \f{s,Ti) is given by (5.19)-(5.20). 

Remark 5.8. Let us point out again that this method can be appHed to 
any LIBOR model driven by a general semimartingalc; indeed, the prop- 
erties of Levy processes are not essential in the proofs of Lemma 5.4 or 
Proposition 5.7. If we start with a LIBOR model driven by a general semi- 
martingale then the structure of this semimartingalc will be "transferred" to 
the first variation process, and hence to the dynamics of the strong Taylor 
approximation. 

6. Numerical illustration 

The aim of this section is to demonstrate the accuracy and efficiency of the 
Taylor approximation scheme for the valuation of options in the Levy LIBOR 
model compared to the "frozen drift" approximation. We will consider the 
pricing of caps and swaptions, although other interest rate derivatives can 
be considered in this framework. 

We revisit the numerical example in Kluge (2005, pp. 76-83). That is, we 
consider a tenor structure Tq = 0, Ti = ^, T2 = 1 . . . , Tiq = 5 = T*, constant 
volatilities 



X{s 


Ti) 


= 0.20 


A(s,T2) 


= 0.19 


X{s 


T3) 


= 0.18 


X{s 




= 0.17 




= 0.16 


X{s 


n) 


= 0.15 


X{s 


Tr) 


= 0.14 


X{s,Ts) 


= 0.13 


X{s 


T9) 


= 0.12 



and the discount factors (zero coupon bond prices) as quoted on February 
19, 2002; cf. Table 6.1. The tenor length is constant and denoted by 5 = |. 

The driving Levy process H is a normal inverse Gaussian (NIG) process 
with parameters a = S = 1.5 and /x = /3 = 0. We denote by fi^ the random 
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measure of jumps of H and by z/(dt,dx) = F{Ax)(\t the Pr* -compensator 
of ix^ , where F is the Levy measure of the NIG process. The necessary 
conditions are satisfied because M = a, hence |A(s,rj)| = 1.44 < a 

and A(s, Tj) < f , for all i G {1, . . . , 9}. 

The NIG Levy process is a pure-jump Levy process and, for = 0, has 
the canonical decomposition 

H = j j x{fi" -iy)ids,dx). (6.1) 

R 

The characteristic function of the NIG distribution is 

'PHiiu) = exp ^6a — 5\/ oP' + (6.2) 
for all II € M; moreover, from (2.6), (2.7) and (6.2), we get that 
k(u) = y"(e"^ - 1 - ux)F{dx) = log]E[exp(tiiJi)] 

= Sa- Wa'^ - u'^, (6.3) 
for all u G C with \^u\ < a. 

6.1. Caplets. The price of a caplet with strike K maturing at time Tj, using 
the relationship between the terminal and the forward measures cf. (3.6), 

can be expressed as 

Co{K, Ti) = SB{0, Ti+i) Ep^^^^ [{L{Ti, Ti) - K)+] 

rdPr 

= SBiO, r,+i) Ep,^ 1^^^ im, Ti) - K)- 

N 

= (55(0, r,) Ep^^ [ n (1 + '^^(^^' ^')) (^(^^' - ^)^] • (6-^) 

l=i+l 

This equation will provide the actual prices of caplets corresponding to 
simulating the full SDE for the LIBOR rates. In order to calculate the first- 
order Taylor approximation prices for a caplet we have to replace L{-,T.) in 
(6.4) with L^(-,r.). Similarly, for the frozen drift approximation prices we 
must use L°(-,T.) instead of L{-,T.). 

The main difficulty when simulating the full SDE for the dynamics of LI- 
BOR rates is that each LIBOR rate depends on all previous LIBOR rates, 
see equations (3.7)-(3.10) or (4.6)-(4.8). Hence, we are facing a system of 
highly nested SDEs that need to be simulated simultaneously. Indeed, the 
"dependence" structure of LIBOR rates takes the form of a (lower) triangu- 
lar matrix shown in Table 6.2. 

The strong Taylor approximation method navigates around this problem, 
since each approximate LIBOR rate depends on all the Taylor approxima- 
tions of the previous LIBOR rates, which arc independent of each other. 
Hence, the resulting SDEs are not nested and can be simulated in parallel, 
see equations (5.18)-(5.20) and (5.8)-(5.9). In other words, we can simulate 
L'^{-,Ti), without having to simulate all previous N — {i + 1) approximate 
LIBOR rates. This also means that the rate of convergence of numerical 



TAYLOR APPROXIMATION OF SDE AND LEVY LIBOR MODEL 



17 







L{t,TN-2) .. 


. L{t,TN-2) 


L{t,TM-i) .. 


. L{t,TN-i) L{t,TN-i) 


... L{t,TN) .. 


. L{t,TN) L{t,TN) L{t,TN) 


... L{t,Ti) .. 


. L{t,TN-3) L{t,TN-2) L{t,TN-i) L{t,TN) 


Table 6.2. 


Matrix of dependencies for LIBOR rates 



schemes for U{-.Tj) is usually higher than for the original system. The nu- 
merical experiments show that the Taylor approximation method reduces 
the computational complexity of the model, while maintaining the accuracy 
of prices. 

We will compare the performance of the strong Taylor approximation 
relative to the frozen drift approximation in two ways: (1) we compute the 
prices of caplets by simulating separately the full SDE, the strong Taylor 
approximation SDE and the frozen drift SDE and compare the prices in 
terms of their implied volatilities; (2) in order to get a better picture of the 
performance of the two methods, we simulate the difference between the 
full SDE and the strong Taylor or the frozen drift prices; this resembles a 
variance reduction technique in Monte Carlo methods and results in a very 
low standard error of the Monte Carlo simulations. 

The implied volatility surface corresponding to the full SDE prices can be 
seen in Table A.l. In Figure 6.1 we present the difference in implied volatil- 
ity between the full SDE prices and the frozen drift prices, while in Figure 
6.2 we show the difference in implied volatility between the full SDE prices 
and the strong Taylor prices. One can immediately observe that the strong 
Taylor approximation method performs much better than the frozen drift 
approximation; the difference in implied volatilities is very low across all 
strikes and maturities. Indeed, the difference in implied volatility between 
the full SDE and the strong Taylor prices lies always below the 1% threshold, 
which deems this approximation accurate enough for practical implementa- 
tions. On the contrary, the difference in implied volatilities for the frozen 
drift approximation exceeds the 1% level for in-the-money options. 

The implied volatilities for all three methods and the standard error of the 
Monte Carlo simulation for the full SDE prices are reported in Appendix A. 
Note that we have used 1 million paths and 200 time steps in all the Monte 
Carlo simulations. 

However, the standard error in these Monte Carlo simulations is relatively 
high, cf. Table B.4. In order to reduce this error and get a more accurate 
picture of the performance of the methods we simulate the differences be- 
tween caplet prices in the different methods and get numerical results with 
relatively low Monte Carlo error. Figure 6.3 exhibits the difference between 
the strong Taylor approximation caplet prices and the full SDE prices, while 
Figure 6.4 exhibits the difference between the frozen drift and the full SDE 
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Figure 6.1. Difference in implied volatilities between the 
full SDE and the frozen drift prices. 



0.15^ 




Figure 6.2. Difference in implied volatilities between the 
full SDE and the strong Taylor prices. 



caplet prices. The respective standard errors for the Monte Carlo simula- 
tion are reported in Tables A. 5 and A. 6. We can immediately observe that 
the difference between the strong Taylor and the full SDE prices is signifi- 
cantly smaller than the difference between the frozen drift and the full SDE 
prices. Moreover, the respective standard error is several orders of magni- 
tude smaller, pointing again to the superior performance of the strong Taylor 
approximation relative to the frozen drift approximation. 

6.2. Swaptions. Next, we will consider the pricing of swaptions. Recall 
that a payer (resp. receiver) swaption can be viewed as a put (resp. call) 
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0.07 



0.02 



Figure 6.3. Difference in caplet prices between tlie strong 
Taylor and full SDE method. 




0.07 



0.02 



Figure 6.4. Difference in caplet prices between the frozen 
drift and full SDE method. 



option on a coupon bond with exercise price 1; cf. section 16.2.3 and 16.3.2 
in Musiela and Rutkowski (1997). Consider a payer swaption with strike 
rate K, where the underlying swap starts at tini6 Ti and. matures at 
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{i < m < N). The time-Tj value is 



ST,(if,7i,r^)= 1- J2 ^kB{Ti,Tk) 



k=i+l 



m / k—\ ^ 

1- E «niT7i^ 




where 



Cfc 



K, i + l<A;<m-l, 
1 + X, k = m. 



(6.5) 



(6.6) 



Then, the time-0 value of the swaption is obtained by taking the Pj;- 
expectation of its time-Tj value, hence 

So = MK,Ti,Tm) 
= B{0,Ti) Ep^^ 
= B{0,T^) 

xEp,, 



N / ™ / 1 \ 

n(i.«(..,.,))(.-Ej=.nrT«p^) 



= B(o,r*)]Ep, 



k='i 



N 



l=k 



(6.7) 



where Ci := —1. Once again, this equation will provide the actual prices of 
swaptions corresponding to simulating the full SDE for the LIBOR rates. 
In order to calculate the first-order Taylor approximation prices we have to 
replace L(-,T.) with L'^(-,r.), and for the frozen drift approximation prices 
we must use L^{-,T.) instead of L{-,T.). 

We will price eight swaptions in our tenor structure; we consider 1 year 
and 2 years as option maturities, and then use 12, 18, 24 and 30 months as 
swap maturities for each option. Similarly to the simulations we performed 
for caplets, we first simulate the prices of swaptions using all three methods, 
and then simulate the differences in prices of swaptions between the full SDE 
and the strong Taylor method, and the full SDE and the frozen drift method, 
in order to reduce the standard error of the Monte Carlo simulation. 

The prices of swaptions corresponding to the simulation of the full SDE for 
the LIBOR rates are reported in Table B.l, while the prices corresponding to 
the strong Taylor method and the frozen drift method arc reported in Tables 
B.2 and B.3 respectively. Once again we observe that the strong Taylor 
method is performing very well across all strikes, option maturities and 
swap maturities, while the performance of the frozen drift method is poor 
for in-the-money swaptions and seems to be deteriorating for longer swap 
maturities. This observation is in accordance with the common knowledge 
that the frozen drift approximation is performing worse and worse for longer 
maturities. 
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Figure 6.5. Difference in swaption prices between tlie 
strong Taylor and full SDE method. 



0.4 




Figure 6.6. Difference in swaption prices between the 
frozen drift and full SDE method. 

Next, we have simulated the differences in swaption prices between the 
full SDE and the strong Taylor method and the full SDE and the frozen 
drift method; the differences in swaption prices are shown in Figures 6.5 
and 6.6 respectively. The standard errors of the corresponding Monte Carlo 
simulations are reported in Tables B.5 and B.6 respectively. The performance 
of the strong Taylor method is clearly superior in comparison to the frozen 
drift method, while we can observe once again that the standard error of the 
full SDE-strong Taylor simulation is significantly smaller than the standard 
error of the full SDE-frozen drift simulation. 
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Appendix A. Caplet prices, standard errors and implied Black 

VOLATILITIES. 

In this appendix we report the impHed Black volatiHties and the standard 
errors of the Monte Carlo simulation for the caplet prices corresponding to 
the full SDE for the LIBOR rates. Wc also report the implied volatilities of 
the strong Taylor approximation and the frozen drift approximation caplet 
prices. The standard errors in these cases are omitted for the sake of brevity; 
we point out though that the standard errors are very similar in all three 
methods. 



2.5% 3.0% 3.5% 4.0% 4.5% 5.0% 5.5% 6.0% 6.5% 7.0% 



T7 



27.27 22.11 19.09 18.88 20.93 23.35 25.59 27.59 29.36 30.88 

26.09 21.71 19.53 18.49 18.36 18.87 19.73 20.71 21.67 22.59 

21.90 19.31 18.12 17.62 17.58 17.88 18.35 18.92 19.51 20.09 

22.22 19.07 17.65 17.01 16.77 16.77 16.94 17.21 17.55 17.91 

19.84 17.42 16.42 15.99 15.82 15.83 15.94 16.14 16.38 16.64 

18.82 16.47 15.51 15.09 14.90 14.85 14.89 15.01 15.17 15.34 

17.24 15.22 14.41 14.07 13.91 13.88 13.92 14.01 14.14 14.28 

15.75 14.03 13.36 13.07 12.93 12.88 12.90 12.95 13.05 13.15 

11.90 12.27 12.10 11.96 11.88 11.85 11.87 11.93 12.01 12.10 



Table A.l. Implied Black volatilities for full SDE prices. 



2.5% 3.0% 3.5% 4.0% 4.5% 5.0% 5.5% 6.0% 6.5% 7.0% 



27.23 22.10 19.09 18.88 20.93 23.35 25.59 27.59 29.36 30.88 

26.02 21.68 19.52 18.48 18.36 18.87 19.73 20.71 21.67 22.59 

21.82 19.28 18.11 17.61 17.58 17.87 18.35 18.92 19.51 20.09 

22.21 19.06 17.65 17.01 16.77 16.77 16.94 17.21 17.55 17.91 

19.93 17.46 16.44 16.00 15.83 15.83 15.95 16.14 16.38 16.64 

19.00 16.56 15.54 15.10 14.91 14.86 14.90 15.01 15.17 15.34 

17.41 15.29 14.44 14.08 13.92 13.88 13.92 14.01 14.14 14.28 

15.84 14.07 13.38 13.08 12.94 12.89 12.90 12.96 13.05 13.15 

11.90 12.27 12.10 11.96 11.88 11.85 11.87 11.93 12.01 12.10 



Table A. 2. Implied Black volatilities for strong Taylor prices. 



2.5% 3.0% 3.5% 4.0% 4.5% 5.0% 5.5% 6.0% 6.5% 7.0% 



25.74 
24.21 
20.54 
20.45 
17.65 
16.13 
15.32 
14.45 
11.90 



21.76 

21.04 
18.89 
18.33 
16.64 
15.53 
14.54 
13.59 
12.27 



19.00 

19.31 
17.97 
17.36 
16.12 
15.15 
14.15 
13.21 
12.10 



18.85 
18.40 
17.55 
16.87 
15.84 
14.91 
13.94 
13.00 
11.96 



20.91 

18.31 
17.54 
16.69 
15.73 
14.80 
13.84 
12.89 
11.88 



23.34 
18.84 
17.85 
16.72 
15.77 
14.78 
13.83 
12.86 
11.85 



25.58 
19.71 
18.33 
16.90 
15.90 
14.85 
13.89 
12.88 
11.87 



27.59 
20.69 
18.90 
17.19 
16.10 
14.97 
13.98 
12.94 
11.93 



29.35 
21.66 
19.50 
17.52 
16.35 
15.14 
14.11 
13.03 
12.01 



30.88 
22.58 
20.08 
17.89 
16.62 
15.32 
14.26 
13.14 
12.10 



Table A. 3. Implied Black volatilities for frozen drift prices. 



TAYLOR APPROXIMATION OF SDE AND LEVY LIBOR MODEL 



23 



2.5% 3.0% 3.5% 4.0% 4.5% 5.0% 5.5% 6.0% 6.5% 7.0% 

0.288 0.274 0.240 0.175 0.116 0.080 0.057 0.042 0.032 0.025 

0.448 0.436 0.412 0.366 0.299 0.230 0.175 0.135 0.107 0.087 

0.502 0.488 0.458 0.406 0.338 0.268 0.209 0.163 0.129 0.103 

0.587 0.575 0.553 0.514 0.456 0.386 0.318 0.257 0.206 0.166 

0.594 0.583 0.560 0.520 0.461 0.393 0.326 0.265 0.213 0.172 

0.619 0.610 0.592 0.558 0.507 0.444 0.376 0.312 0.256 0.208 

0.600 0.592 0.575 0.542 0.493 0.432 0.366 0.304 0.248 0.202 

0.592 0.587 0.575 0.549 0.506 0.450 0.387 0.324 0.267 0.218 

0.554 0.551 0.540 0.517 0.477 0.423 0.361 0.301 0.245 0.198 

Table A. 4. Standard error xlO~^ of the Monte Carlo sim- 
ulation for the full SDE prices. 





2.5% 


3.0% 


3.5% 


4.0% 


4.5% 


5.0% 


5.5% 


6.0% 


6.5% 


7.0% 




0.006 


0.019 


0.047 


0.080 


0.057 


0.038 


0.027 


0.019 


0.015 


0.012 




0.011 


0.031 


0.070 


0.120 


0.211 


0.193 


0.156 


0.121 


0.094 


0.075 


T3 


0.023 


0.060 


0.120 


0.182 


0.316 


0.299 


0.256 


0.209 


0.168 


0.135 


Ti 


0.004 


0.010 


0.020 


0.033 


0.043 


0.062 


0.060 


0.054 


0.046 


0.039 


T5 


0.031 


0.078 


0.152 


0.237 


0.304 


0.533 


0.513 


0.463 


0.403 


0.343 




0.048 


0.123 


0.245 


0.394 


0.529 


0.611 


0.982 


0.925 


0.832 


0.726 




0.047 


0.118 


0.234 


0.375 


0.500 


0.576 


0.912 


0.861 


0.776 


0.679 




0.015 


0.039 


0.081 


0.135 


0.188 


0.226 


0.366 


0.355 


0.326 


0.290 



Table A. 5. Standard error xlO ^ of the Monte Carlo sim- 
ulation for the difference in caplet prices between the strong 
Taylor and the full SDE method. 





2.5% 


3.0% 


3.5% 


4.0% 


4.5% 


5.0% 


5.5% 


6.0% 


6.5% 


7.0% 




0.002 


0.006 


0.014 


0.024 


0.017 


0.012 


0.008 


0.006 


0.004 


0.004 




0.003 


0.007 


0.016 


0.027 


0.047 


0.043 


0.035 


0.027 


0.021 


0.017 


n 


0.004 


0.009 


0.019 


0.028 


0.050 


0.047 


0.040 


0.033 


0.027 


0.021 




0.004 


0.010 


0.021 


0.034 


0.045 


0.078 


0.075 


0.066 


0.057 


0.048 


T5 


0.006 


0.015 


0.028 


0.044 


0.057 


0.098 


0.094 


0.085 


0.074 


0.063 




0.005 


0.013 


0.025 


0.041 


0.055 


0.063 


0.101 


0.095 


0.086 


0.075 


T7 


0.004 


0.010 


0.020 


0.032 


0.043 


0.049 


0.078 


0.074 


0.067 


0.058 


Ts 


0.002 


0.004 


0.009 


0.014 


0.020 


0.024 


0.039 


0.038 


0.035 


0.031 



Table A. 6. Standard error xlO^'^ of the Monte Carlo sim- 
ulation for the difference in caplet prices between the frozen 
drift and the full SDE method. 



Appendix B. Swaption prices and standard errors. 

In this appendix we report the swaption prices for the three methods and 
the standard error of the Monte Carlo simulation for the swaption prices 
corresponding to the full SDE for the LIBOR rates; we point out that once 
again the standard errors are very similar in all three methods. We also 
report the standard errors corresponding to the simulation of the difference 
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in swaption prices. In all the simulations we have used 1 million paths and 
200 time steps. 
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Table B.l. Swaption prices for the full SDE method. 
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Table B.2. Swaption prices for the strong Taylor method 
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Table B.3. Swaption prices for the frozen drift method. 
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Table B.4. Standard error of the Monte Carlo simulation 
for the full SDE prices. 
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Table B.5. Standard error xlO~ of the Monte Carlo simu- 
lation for the difference in swaption prices between the strong 
Taylor and the full SDE method. 
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Table B.6. Standard error x 10-2 of the Monte Carlo simu- 
lation for the difference in swaption prices between the frozen 
drift and the full SDE method. 
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